Spin-Orbit Qubits of Rare-Earth-Metal Ions in Axially Symmetric Crystal Fields 
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Contrary to the well known spin qubits, rare-earth qubits are characterized by a strong influence 
of crystal field due to large spin-orbit coupling. At low temperature and in the presence of resonance 
microwaves, it is the magnetic moment of the crystal-field ground-state which nutates (for several 
/j,a) and the Rabi frequency Qr is anisotropic. Here, we present a study of the variations of Qr(Ho) 
with the magnitude and direction of the static magnetic field Ho for the odd 167 Er isotope in a 
single crystal CaW04:Er 3+ . The hyperfine interactions split the Qr(Hq) curve into eight different 
curves which are fitted numerically and described analytically. These "spin-orbit qubits" should 
allow detailed studies of decoherence mechanisms which become relevant at high temperature and 
open new ways for qubit addressing using properly oriented magnetic fields. 



PACS numbers: 71.70.Ej,75.10.Dg,76.30.Kg, 03.67.-a 

Worldwide studies on possible implementation of a 
spin-based quantum computer showed the significance 
of scalable single electron spin-qubits in e.g. quantum 
dots [1, 2, 3]. Are also interesting atomic and nuclear 
spins 1/2 dispersed in solid state matrices with coherence 
times persisting at relatively high-temperature [4, 5, 6] 
or a new class of potentially scalable systems which re- 
cently emerged with the observation of Rabi oscillations 
[7] in the single molecular magnet V15 [8], followed by 
electronic-spin coherence studies on Fes [9] and Fe4 [10]. 
Spin 1/2 qubits are isotropic or quasi- isotropic and their 
Rabi frequency flf> = .gyUs/i m ,„/2/i depends only on the 
amplitude of the linearly polarized microwave field h mw 
(h is the Planck constant, /U_b is the Bohr magneton and 
9~2). 

With rare-earths ions in crystals [11], the situation is 
very different. Contrary to qubits with magnetic mo- 
ments simply proportional to the spin (single electrons, 
3d transition metal ions with zero or quenched orbital 
moment), the orbital magnetic moment of a rare earth 
ion is large and not reduced by the crystal-field, this last 
being much weaker than the spin-orbit coupling. In gen- 
eral case this leads to a strong magnetic anisotropy and, 
as it will be shown in this paper (see preliminary results 
in [11]), to a strong dependence of the Rabi frequency on 
the directions of the microwave and static applied fields. 
The amplitude of the total angular moment J = L + S 
is generally much larger than 1/2 leading to magnetic 
moments as large as 10/is (Ho 3+ ) allowing spin manipu- 
lations in low driving fields. For odd isotopes, hyperfine 
interactions of 4/ electrons with the rare earth nuclear 
spin can be large enough so that nuclear spins are entan- 
gled with the total angular moment [12] and modify the 
Rabi frequency [8]. 

In this paper we concentrate on the anisotropy of the 
Rabi frequency of a single crystal CaW04:Er 3+ , a system 



which presents an in-plane magnetic anisotropy favoring 
the observation of Rabi oscillations and several Erbium 
isotopes (odd and even) allowing the study of coherent 
dynamics with and without nuclear spin. In the first 
part we give a theoretical description (numerical and 
analytical) of what we may call "spin-orbit qubits" and 
in the second one we present the experimental results 
confirming quantitatively the theoretical predictions and 
allowing to understand the effect of strong spin-orbit 
coupling on Rabi oscillations. 

The crystals of CaWO 4 :0.05% Er 3+ used in this study, 
obtained by the Chokralsky method, are characterized 
by a body centered tetragonal scheelite-type structure 
(I Ai/ a space group) with the lattice constants a = b = 
0.524 nm and c = 1.137 nm [13]. The Er 3+ ions sub- 
stitute for the Ca 2+ ions inside the eight-fold oxygen 
surroundings, charge excess being compensated by sup- 
plementary substitutions of Na + ions in the crystal. 
The Hamiltonian operating in the space of states of the 
ground 4 /i5/2 multiplct of a single Er 3+ ion in CaW04 
(S4 point symmetry) is given by: 

H = H cf + H hf +H z (1) 

where H c f is the Hamiltonian of a single rare-earth ion, 
reduced by the host-matrix crystal-field and 4/ symme- 
try: 

H cf = a.jB 2 O a 2 +l3j(B 4 O° 4 +BtOi + B^O^) 

+lJ {BlOl + BtOl + B^O^) (2) 

The O™ are the Stevens equivalent operators, ctj, 
(3 j, 77 the Stevens coefficients [14] and the are the 
crystal-field parameters. For CaW0 4 :Er 3+ , £^=231, 
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cm 1 in the Cartesian system of coordinates 
with the z-axis parallel to the c-axis. The second 
term in (1), Hhf = AjJ.I, is the magnetic hyperfme 
interaction for the 167 Er isotope (I = 7/2, 22.9% 
abundance, hyperfme constant Aj = -125MHz [15]), 
quadrupolar couplings are negligible. The third term 
in (1), Hz = gjHsJ-Ho (gj = 6/5 is the Lande factor 
of Er 3+ ), is the static Zeeman interaction. We neglect 
here dipole-dipole interactions between the highly 
diluted Er 3+ ions and super-hypcrfine interactions 
with nuclear magnetic moments of the 183 W isotope. 
These interactions induce negligible changes of Rabi 
frequencies, but may be principal sources of decoherence 
[16, 17, 18, 19]. The crystal-field of CaW0 4 leads to 
strong in-plane anisotropy of Er 3+ (a — b plane) [15] 
and effective g-factors g^ = g c =1.2, g± = g a = g b = 8.4 
[20]. The Rabi frequency can be calculated either in 
the laboratory frame (LF) by computing the amplitude 
of probability of transition induced by the microwave 
field h mw perpendicular to the static field Hq or in 
the rotating frame (RF) by computing the splitting at 
avoided level crossing. 

In the LF, we start with the Hamiltonian (1) of 
a single 167 Er 3+ ion. By diagonalization of the 128 
dimension Hilbert space electro-nuclear Hamiltonian 
(J=15/2cS>I=7/2), we calculate the Rabi frequencies as 
a function of the direction of the static field Hq us- 
ing the expression flu = \ gj^BJpk-h m w | /2h where 
Jpk = (4>p{Ho)\J\4>k(Ho)) is a vector varying with the 
static field Hq and (j)p,k(Ho) are the wave functions re- 
sulting from the diagonalization of (1). In general, the 
anisotropy of flu depends on both the directions of Hq 
and h mw . These expressions for spin-orbit qubits in 
a crystal-field generalize the one of simple spin qubits, 
VLr = g a HBhf nw /2h [21] (g a is the effective g-factor in 
the direction a of the electromagnetic wave polarization 
h a ) 

In the RF, we use the time-dependent effective Hamil- 
tonian of the ground states doublet : 



H(t) = f i B H Q [g]S + I[A}S + f i B h mw [g}Scos{2irft) (3) 

/ is the microwave frequency, [g] the g-tensor, [A] = 
[g]Aj/gj the hyperfme tensor. The NMR paradigm ap- 
plied to EPR [22, 23] allows to transform (3) into a 
time independent effective Hamiltonian. The Rabi fre- 
quencies are then given by avoided levels splitting. The 
Rabi frequencies of the 1=7/2 167 Er 3+ isotope are derived 
from the RF energy spectrum obtained numerically for 
,ffo||a = b, Fig.l. The effective electro-nuclear spin-states 
[1/2, mi) and | — 1/2, mj) coupled by the electromagnetic 
wave ("dressed states") form several avoided level cross- 
ings of which splitting is hfln (decreasing rapidly when 
the microwave field rotates from the easy a — b plane to 
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FIG. 1: (color online) Energy spectrum of 167 Er 3+ with 7=7/2 
calculated numerically in the rotating frame. The detuning 
field is, by definition, the difference between the applied field 
Ho and the resonant field of the 7=0 isotope. The Rabi fre- 
quencies are equal to the splittings at avoided level crossings 
(dressed states). Allowed EPR transitions (Ami = 0) are 
gapped and their anisotropy for different directions of the mi- 
crowave polarization is shown in the inset. 



the hard c-axis). In the case of the 1 = isotope, when 
Hq makes an angle <f> with c-axis in the x-z plane and 
h m w is in the same plane, the result can be analytical 
and given by: 



9±9\\ 
'2hg(ct>) 



(4) 



where g(4>) = We^ cos 2 <fi + g\ sin 2 </>. Much more ex- 
tended analytical results will be given in a forthcoming 
publication [24] showing that both LF and RF method 
give, in all cases, quite the same results. In the following 
we will use the LF method which is more convenient for 
precise calculations. 

The EPR measurements have been performed using a 
conventional Bruker EPR spectrometer E580. This spec- 
trometer works at both, continuous wave (CW) or time 
resolved (TR) mode in the X band. All experiments were 
carried out at the frequency of / = 9.7 GHz. The detec- 
tion of the signal was performed using a cavity work- 
ing with the TEon mode with Q-factor of about 4000 
in CW mode (perfect coupled cavity) and 200 in TR 
mode (over-coupled cavity). The sample-cavity ensem- 
ble is surrounded by a 4 He flux cryostat and the tem- 
perature was controlled by an Oxford instrument ITC 
form 2.5 to 300 K. Two pulse sequences have been used 
i) the usual "spin-echo sequence" which strongly reduces 
the decoherence associated with inhomogeneous CW line- 
width, and ii) the Rabi oscillations technique consisting 
in a spin-echo sequence preceded by an excitation pulse. 
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FIG. 2: (color online) (a) Rabi oscillation recorded on the 
nuclear- spin free Er 3+ isotope in CaW0 4 (H = 325. lmT, 
hmw = 0.13mT) and (b) its Fourier transform, (c) Con- 
tour plot of the Rabi frequency distribution measured in the 
I = 7/2 isotope while Ho tilted by 12° from the c-axis is 
swept from 220mT to 440mT (colored iso-magnetization scale 
increases from blue to red). The eight Rabi frequencies with 
smaller intensity correspond to the transitions (1/2, mi) — > 
(—1/2, mi) of the 1—7/2 isotope (Ami = 0) and the single one 
with larger intensity to the 1—0 isotope. Calculated Qr are 
represented by symbols going from black square (mi = — 7/2 
) to green circle ( mi = 7/2). 



This technique allows to obtain the time evolution of the 
averaged magnetization proportional to (S z (t)), i.e. the 
Rabi oscillations if the timescale is such that coherence 
wins over relaxation (here the z— axis ||#o)- 

Fig. 2(a) and (b) give an example of measured Rabi os- 
cillation and its Fourier transform whereas Fig. 2(c) shows 
the eight electro- nuclear transitions of the 1 = 7/2 iso- 
tope, plus the single one of the nuclear-spin free isotope. 
In this experiments Ho was tilted by <fi = 12° from the 
c-axis in order to resolve the eight electro-nuclear tran- 
sitions . The symbols correspond to the Qr= 6.61, 7.30, 
7.95, 8.57, 9.14, 9.65, 10.1, 10.5MHz calculated from di- 
agonalization of the electro-nuclear Hamiltonian (1) (LF 
method) using the published crystal-field parameters and 
hyperfine constant only [15]. The RF expression (4) gen- 
eralized to the case / ^0 and linear in Aj give the values 
Q R = 6.63, 7.22, 7.82, 8.42, 9.02, 9.6, 10.21, 10.80MHz 
[24]. The small differences between the two methods are 
caused by the neglect of higher terms in Aj. Both sets 
of values agree very well with experimental data. 

Fig. 3 shows the evolution of the Rabi frequency of 
each electro-nuclear transition when the crystal rotates 
from 9 =0° (i?o||c') to 9 =90° (H \\b) (see inset). Above 
9 « 30° (not shown) all the Rabi frequencies ^Ir(itij,9) 
tend asymptotically to the value Qr ~ 15MHz predicted 
by the usual equation Qr = g a fj,Rh^ lw /2h [21] inde- 
pendently of mj. At 9 < 30° these quasi-degenerated 



N 

X 



a 



16 



14 



12 



10 



1 1 1 I 1 1 


' •••• jpp a 

■ a 

■~3W V 


1 =o.it\ 

h =0.26mT 

^ mw 

mw 








1 i 1 i 





10 20 
Angle e (Deg) 



30 



FIG. 3: (color online) Evolution of S}r when the static field 
Ho rotates by an angle 6 around the microwave field hmw 
Colored continuous lines are computed numerically. Inset: 
geometry of the experiment where the sample defined by its 
axes a, b, c is tilted by a small angle (f> ~ 7° around the b=b'- 
axis. The dotted line represents the conventional equation 
Q R = g a fJ , B h mw /2h [21]. 



frequencies VLR(mi,9) split off in a way strongly de- 
pendent on the algebraic value of m/. It is the small 
tilt angle 4> which allows to resolve the transitions with 
different mi (i.e. for <fi =0, Qr does not depend on 
77i/). Non-magnetic isotopes show a similar variation 
laying between the blue and green triangles with A^ = 
Qr(Ho\\c')/Qr(H _L c') « 0.77. In these experiments a 
rotation of the resonance field Ho implies a change of its 
amplitude in order to keep at resonance, and this rises 
the question to know whether the observed variations 
of Qn(mi,9) are mainly due to changes in direction or 
amplitude of the vector field Hq- In order to answer this 
question, we plot in Fig. 4 the evolution of Qr vs j| H ||. 
The large dispersion observed on Qr at given || Ho \\ 
definitively confirms that it is, as expected, a directional 
effect. All the measured Rabi frequencies of Fig. 3 and 
Fig. 4 are very well reproduced by our crystal-field model 
without any fitting parameter using only published val- 
ues of parameters, as described above. In particular, the 
ratio Ar varies from 0.58 to 0.90, whereas, in absence of 
orbital contribution, AR(mi) = 1 for any value of mi. 

Interestingly, if the coordinate system is such 
that the microwave field h mw has only one non- 
zero component (hf nw ), the scalar-product equation 
= l 9j^bJ -hmw | /2h takes the simpler form : 
Q R (Ho,h mw ) = G(a,f3)tx B h^ w /2h where G(a,/3) = 



gj\{(f>i(H )| J a \4>2{Ho ))\ depends on both directions a 
and (3 of the static and microwave fields (a. (3 = 0). 
This expression traduces the fact that, when Ho rotates 
around h mw the crystal-field/Zceman competition varies 
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FIG. 4: (color online) Evolution of Q,r with Qr vs || Ho || 
while the resonant field is rotated. The symbols represent 
the same data and the lines are calculated numerically and 
analytically as in Fig. 3. The dotted line represents the con- 
ventional equation Qr = g a )j,Bh% lw /2h [21]. 



symmetry and temperature-dependent spin-lattice deco- 
herence mechanism; in particular, the transfer of quan- 
tum dynamics from the spin-bath to the nutating sys- 
tem is probably partial implying the existence of resid- 
ual decoherence [16, 17, 18, 19]. (ii) Despite significant 
spin-phonon coupling, relatively long living coherence is 
observed (« 50/xs at 2.5K in CaWC>4 :Er 3+ ) showing 
that the coupling with crystal-field environment is not 
redhibitory, (iii) The magnetic moment generally much 
larger than 1/2 allows spin manipulations in low driving 
field- vectors (amplitude and direction), (iv) RE qubits 
inserted in a semi-conducting film, should be scalable 
with selective addressing (application of weak local field 
pulses created by nano-line current adding algebraically 
to uniform static field) and couplings (controlled carrier 
injection through a gate voltage). 

We thank for financial support the former INTAS-99- 
01839 contract in which the crystals of Er:CaWC>4 were 
synthesized, the CEA, and the European Network of Ex- 
ellence MAGMANET. 



leading rcnormalization of wave functions entailing varia- 
tions of G(a,(3) , i.e. of Hr. As mentioned above, this is 
true only if the crystal is tilted so that h mw is not along a 
crystallographic axis (<fi ^0, inset Fig. 3). In the opposite 
case where h mw is along a crystallographic axis (a-axis 
in Fig. 3) £Ir is isotropic and follows the usual expression 
fin = g a iiBh^ nw j1h [21] represented by a dashed line 
in Fig. 3 and 4. This renormalization of wave functions 
appears clearly in the analytical expressions for the non- 
magnetic and magnetic isotopes: £Ir depends on both g\\ 
and <7j_ if 4> ^0 whereas for <f) =0 it depends on g± only, 
and Qft(mj) associated with the different mj are clearly 
different from each others unless <f> =0 or 7r/2 [24] . 

In conclusion, we have shown with the example of Er 3+ 
ions diluted in the single-crystalline host matrix CaWO.4, 
that crystal-field deeply modifies the coherent quantum 
dynamics of what we call "spin-orbit qubits" . In particu- 
lar, it introduces dramatically large variations of the Rabi 
frequency fin when the static or/and dynamical field (Hq 
or h mw ) deviate from crystallographic axes. This effect, 
directly connected with modifications of the wave func- 
tions of resonant states is amplified or depressed by the 
hypcrfine interaction of the odd isotope (7 =7/2) which 
removes the eight-fold degeneracy of Qr by increasing 
or decreasing it when mi >0 or mi <0. All the ex- 
perimental data, showing a rich pattern of coherent os- 
cillations, are interpreted quantitatively without requir- 
ing any fitting parameter giving a clear picture of this 
new type of crystal- field dependent qubits. These qubits 
are interesting for several reasons (i), the coherent nuta- 
tion of the ground-state magnetic moment deriving from 
crystal-field effects acting on J = L + S (and J + I with 
odd isotopes) is associated with , not yet well studied, 
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